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ABSTRACT 

We compare the efficiency with which 2D and 3D weak lensing mass mapping tech¬ 
niques are able to detect clusters of galaxies using two state-of-the-art mass reconstruc¬ 
tion techniques: MRLens in 2D and GLIMPSE in 3D. We simulate otherwise-empty 
cluster fields for 96 different virial mass-redshift combinations spanning the ranges 
3 x 10 13 h~ 1 M & ^ M V i r ^ 10 15 /i _ 1 Mq and 0.05 ^ z c \ ^ 0.75, and for each generate 
1000 realisations of noisy shear data in 2D and 3D. For each field, we then compute 
the cluster (false) detection rate as the mean number of cluster (false) detections 
per reconstruction over the sample of 1000 reconstructions. We show that both MR¬ 
Lens and GLIMPSE are effective tools for the detection of clusters from weak lensing 
measurements, and provide comparable quality reconstructions at low redshift. At 
high redshift, GLIMPSE reconstructions offer increased sensitivity in the detection 
of clusters, yielding cluster detection rates up to a factor of ~ 10x that seen in 2D 
reconstructions using MRLens. We conclude that 3D mass mapping techniques are 
more efficient for the detection of clusters of galaxies in weak lensing surveys than 2D 
methods, particularly since 3D reconstructions yield unbiased estimators of both the 
mass and redshift of the detected clusters directly. 

Key words: gravitational lensing: weak - cosmology: dark matter - cosmology: large- 
scale structure - galaxies: clusters: general 


1 MOTIVATION 

Current and upco ming large surveys such as the D ark En¬ 
ergy Survey IDE S. iT he Dark Energ y Survey Collaboration! 
120051 1. Euclid (iRefregier et al.l l2010l h and the Large 
Synoptic Survey Telescope (LSST, llvezic et al.l 120081 : 
ILSST Science Collaboration et alj 120091 ') will soon provide 
a large volume of high-quality weak lensing data covering a 
significant fraction of the sky up to redshifts z ~ 2. While 
the primary weak lensing science aims of these projects cen¬ 
tre around the use of two-point statistics to constrain the 
cosmological parameters, there is nonetheless significant in¬ 
terest in reconstructing the mass distribution inferred by 
the lensing shear field. Such maps are obtained essentially 
for free in weak lensing surveys, as mass reconstruction algo¬ 
rithms require no additional data inputs to two-point shear 
statistics. 

In addition to providing a valuable method for visu¬ 
alising the distribution of structure in the Universe, such 
maps can help to facilitate the computation of higher-order 
lensing statistics, in addition to providing a straightforward 
method to detect large structures such as clusters of 
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galaxies. This can provide substantial insights into the 
nature and evolution of nonlinear structure in the Universe, 
and can help to break the degeneracy seen between cosmo¬ 
logical parameters in constraints obtained from two-point 
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Until very recently, weak lensing mass reconstruction efforts 
focused on estimating the two-dimensional projected mass 
density, the convergence k, integrated along the line of sight. 
With the advent of accurate photometric redshift estimation 
methods and high-quality data, it has now become possible 
to con sider reconstructing the fu l l three - dimensional density 
field (ISimon. Taylor fe Hartlapl 120091 : ISimon et al.l 1201 ll : 
ISimonl 120121 ; VanderPlas et al.l 1201 lh . One advantage to 
this approach is that, by directly reconstructing the density 
field, we might expect to be able to directly estimate the 
masses of haloes detected without relying on assumptions 
about the dynamical state of the h alo. In a recent publi¬ 
cation, iLeonard, Lanusse fe Starcfd (I20l4) demonstrated a 
sparsity-b ased approach to 3D l ensing ma ss reconstructions 
(see also ILeonard. Dupe fe Starckl 120129 that is able not 
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only to detect massive clusters of galaxies, but also to 
estimate their masses and redshifts in an unbiased way. 
This represents the first method able to reconstruct the 
density field from weak lensing measurements sufficiently 
accurately to allow direct estimates of the properties of 
dark matter haloes detected. 

Such reconstructions, if sufficiently sensitive, might 
therefore be able to augment and complement optical clus¬ 
ter studies that aim to identify and weigh clusters of galax¬ 
ies by considering the distribution and/o r dynamics of the 
clust e r members (e.g. Girardi_et_alJ 1998|- Lokas&_Mamon 

2003|; iLokas et al.l 120061 : iMamon fc Bouel l2010l ; iFalco et al 


20141 ') . Clusters of galaxies are a very important cosmologi¬ 


cal p robe. Stud yi ng the evolution of the cluster ma ss func¬ 
tion (iRosati, Borgani fe Normanl 120021 ; I Voitl l2005h . allows 
us to constrain both the amplitude of the power spectrum 
at the cluster scale and the linear growth rate of density 
perturbations. In addition, the clustering properties of the 
large-scale distribution of clusters provides direct informa¬ 
tion on the shape and amplitude of th e underl yin g dark 
matter distribution and power spectrum (IBorgani fe Guzzd 
l200ll : iMoscardini, Matarrese fe M 0 II 2 O 0 I] ') . 

In this paper, we focus on two state-of-the-art weak 
lensing mass reconstruction techn iques: GLIMPSE in 3D 
( Dxmnrd^/jm^ 1 4l l and MRLens in 2D 

( Starck, Pires fe RefregieJ 20061 ') The methodology under¬ 
pinning the MRLens algorithm is qualitatively similar to 
that of GLIMPSE: both methods use wavelets, and im¬ 
pose sparsity to regularise the reconstructions and to con¬ 
trol noise peaks. Other reconstruction methods in 2D could, 
of course, be chosen. An obvious example is the aperture 
mass statistic, which is frequently used for 2D mass map¬ 
ping _and_^eak_detectiom H owev er, it was demonstrated 
in iLeonard, Pires fo Starckl (l2012f) that the aperture mass 
statistic is formally identical to the wavelet transform of the 
convergence map, with the added advantages that 


(i) the wavelet transform considers several angular scales 
simultaneously, while application of the aperture mass 
statistic requires the choice of a scale radius, 

(ii) the wavelet transform is b etween 10 to 1 000 times 
faste r depending on the scale (ILeonard, Pires fe Starch 
I2012IL and 


(iii) the wavelet function intrinsically has some desirable 
properties, such as regularity, and compact support, and 
as a compensated filter it does no t require truncation, un- 


like some aperture mass filters (see 

Leonard. Pires & Starckl 

20121; Pires. Leonard & Starck 2012 

, and references therein). 


It was further shown in lPires, Leonard fc Starckl d2012l ~) that 
MRLens produces 2D reconstructions that are more ro¬ 
bust to noise peaks than the aperture mass statistic, and 
are therefore more useful for cosmological studies such as 
those involving peak counts. It was aso shown that MRLens 
outperforms other converge nce reconstruction methods suc h 
as inverse Wiener filtering (IStarck, Pires fc Refregieill2006l h 
For these reasons, and given the similarity to GLIMPSE in 
its wavelet- and sparsity-based approach, MRLens seemed 
the natural choice of algorithm for such a comparison. 

We examine the ability of weak lensing mass reconstruc¬ 
tion techniques in two- and three-dimensions to detect clus¬ 
ters as a function of the cluster mass and redshift. Using 
simulated data, we consider 96 samples in virial mass and 


redshift j = {1..96}, and for each generate 1000 

noise realisations. We generate reconstructions of these noisy 
fields using MRLens and GLIMPSE, and consider both the 
probability of detecting a real cluster of a given mass and 
redshift (the true detection rate), and the probability of find¬ 
ing a peak in a given field that arises due to noise (the false 
detection rate). Comparing the true detection rate in 2D and 
3D reconstructions at the same peak signal-to-noise thresh¬ 
old allows us to probe any intrinsic differences in the lens¬ 
ing signal-to-noise when the full 3D information is retained 
versus when this information is projected along the line of 
sight. However, naively we might expect the false detection 
rate to scale with the number of pixels in the reconstruction, 
which is always larger in 3D than in 2D. The ideal recon¬ 
struction algorithm will maximise the true detection rate 
and minimise the false detection rate, and we examine how 
these compare between the two algorithms, and scale as a 
function of user-specified algorithm parameters. 

This paper is organised as follows: In § [2] we provide 
a brief overview of weak lensing mass reconstructions, and 
outline the MRLens and GLIMPSE algorithms. In § [ 3 ] we 
describe the suite of cluster simulations, and present the 
results of our study as described above. We conclude with 
a brief discussion of our results and their implications for 
cluster and weak lensing surveys in § [4] 


2 MASS RECONSTRUCTION METHODS 


In what follows, we assume that the vector 

7 <‘> - (4*>«<‘>..4>)‘ m 

contains the complex ellipticities of galaxies in a tomo¬ 
graphic weak lensing survey binned on a grid with N g grid 
points. Sources are binned together on the angular grid only 
if they belong to the same redshift bin, which is denoted 
by the superscript ^ in equation © above. We can express 
the full three-dimensional weak lensing information as the 
vector 7 = ( 7 , 7 ^, where N z is the number of 

redshift bins in the survey. 

The aim of mass reconstruction methods is to relate 
the measured shear 7 to the underlying density distribu¬ 
tion. Typically, in 2D, we aim to recover the projected, di¬ 
mensionless surface density (the convergence) F ' 1 ' 1 from the 
tomographic shear 7 ^, while in 3D we seek the dimension¬ 
less 3D matter overdensity <5. We can relate the tomographic 
shear 7 ^ to the convergence k! l ' > of the lens by 

7 w (0) = i J d 2 e'v{e - e')n (i) (e') , (2) 


where 9 = 6 1 +162 represents the angular coordinate on the 
sky, represented in complex notation, 


V{ 6 ) 


1 

WY ’ 


( 3 ) 


and the asterisk * represents complex conjugation. 

Equation © represents a linear mapping between the 
shear and the convergence, and can be expressed in matrix 
notation as 7 ^ = + nSj\ where nty is a vector of 

intrinsic source ellipticities and the linear transform P 7K is, 
for i ^ m: 


> 1 _ 
7«J £m — 


7T [91 


( 4 ) 
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where A represents the solid angle of the angular grid pixels, 
and l and m give the x and y coordinates of the pixel on 
the grid. When £ = m, we set P 7re = 0. 

The underlying convergence, though representing a two- 
dimensional projected density, is actually itself a three- 
dimensional field: at each point r in three-dimensional space, 
we can compute a convergence n(r), by considering the re¬ 
lationship of the lensing convergence to the underlying mat¬ 
ter overdensity S(r). Note that we often parameterise the 
3D position vector r in observations by an angular position 
9 and a redshift 2 . This matter overdensity, or density con¬ 
trast, is defined as S(9,z) = ( p{9,z ) — p(z))/p(z), where 
p(9, z ) is the density at angular position 9 and redshift 2 
and p( 2 ) is the mean matter density at redshift 2 , and is 
related to the convergence by 

w (g ,*)= 3 ^' n 9 M r dw ,fKW)fK(*-w')S[fKW)Q,w'] ' (5) 

2c 2 Jo Ik{w) a{w') 

where Ho is the hubble parameter, Qm is the matter density 
parameter, c is the speed of light, a(w) is the scale parameter 
evaluated at comoving distance w, and 

( K~ 1/2 sin (K 1/2 w), K > 0 

Jk O) = < w, K = 0 , (6) 

[ (—-ff) _ 1 / 2 sinh([—it'] 1 / 2 ui) K < 0 


gives the comoving angular diameter distance as a func¬ 
tion of the comoving distance and the curvature, K , of 
the Universe. fx(w)9 in the argument of 5 in equation lf5|l 
therefore gives the transverse comoving distance. If we now 
consider that we have a tomographic measurement of the 
convergence, k. = (k/ 1 ^ , k) 2 \ ..., k^ N z ^), obtained from the 
tomographic shear 7 , the above relationship can also be 
expressed in matrix notation as ft = Q 5, where Q repre¬ 
sents the line of sight convolution operation in equation (|5|). 
and is given in f ull in Simon. Taylor fe Hartlapl (120091 ) and 
ILeonard. Lanusse fe Starckl ( 2014 1. and 6 now represents a 
binned map of the density contrasijj. For a given set of shear 
measurements 7 , we therefore can write 


7 = P 7 «re + n 7 , 

— P 7K Q<5 4" rt 7 , 


( 7 ) 

(8) 


where, following the conventions used previously, ra 7 = 
(nly \ n 7 2 \ ..., n^y Nz ^) represents the error in the shear mea¬ 
surements. While there are various sources of error in shear 
measurements, for the purposes of this paper we consider 
only noise arising from the intrinsic shapes of galaxies, which 
is typically taken to be a Gaussian distribution of width 
07 ~ 0.2 - 0.3. 

The task of mass reconstruction methods is the inver¬ 
sion of either equation © in 2D or equation © in 3D. 
Note that the formalism presented above assumes the ab¬ 
sence of intrinsic alignments, and this assumption is prop¬ 
agated throughout the analysis. A full assessment of the 
impact of intrinsic alignments on the reconstruction quality 
is the subject of future work. 

We outline below the methods underlying two 
state-of-the-art mass reconstruction techniques: MRLens 


1 Note that the redshift binning in S does not have to be the 
same as that of the convergence and shear; i.e. the matrix Q can 
be rectangular rather than square. 


( Starck J _Pires_& i ^ 1 efregie^_ 2006|) in 2D and GLIMPSE 
I Leonard. Lanusse fc Starckl 2014 1 in 3D. 


2.1 2D mapping with MRLens 

In 2D, our aim is to solve Equation © for each tomo¬ 
graphic redshift bin and estimate the convergence from 
noisy shear measurements. In the Fourier plane, the con- 
vergence can be compu ted from the shear components as 
jKaiser fo Sauireslll993l ) 

/« w (k) = Pi(k) 7 1 w +P 2 (k) 7 2 w , (9) 


where the hat symbol denotes Fourier transforms and we 
have defined k 2 = k 2 + k% and 


Pi(k) = 


k\ - k\ 
k 2 


P 2 (k) 


2fcifc 2 

k 2 


( 10 ) 


with Pi(k) = 0 when k 2 = k 2 , and P 2 (k) = 0 when fci = 0 
or fc 2 = 0. These conditions on Pi and P 2 when fci = fc 2 = 0 
correspond to an indetermination of the mean value of the 
convergence field re, known as the mass-sheet degeneracy. 
Applying an inverse Fourier transform on re estimated by 
equation © yields an estimate of the convergence. As this 
estimate is typically very noisy, some type of filtering or reg- 
ularisation is usally included in a reconstruction algorithm 
in order to minimise the contaminating effect of noise on the 
resulting reconstruction. Gaussian and Wiener filtering are 
common choices, but may not be optimal, particularly since 
Wiener filtering is only optimal when both the signal and 
noise are Gaussian-distributed, and their power spectra are 
known. 

MRLens is a Multi-Resolution Entropy filtering method 
using a wavelet based prior for the entropy to regularise the 
solution. For full details of the MRLens a l gorith m, the reader 
is referred to IStarck, Pires fe Refregied (I2006h : briefly, this 
method assumes that the reconstructed convergence can be 
represented sparsely when projected onto a well-chosen dic¬ 
tionary or set of functions (basis or frame). By this we mean 
that most of the information content can be captured by 
a small number of coefficients in the chosen dictionary or 
frame. For example, a periodic signal would be sparse in the 
Fourier domain, as it would be completely represented by a 
small number of Fourier coefficients. In the present case, we 
choose wavelets. 

In the case of Gaussian noise, one can construct a mul¬ 
tiresolution support in the wavelet domain by identifying 
as “significant” those wavelet coefficients 117 ,r,m that verify 
|«h,q m | > koj, where aj is the noise standard deviation at 
a particular wavelet scale j, and the indices [£, m] again de¬ 
note the pixel location in the 2D convergence map@. The 
noise standard deviation can be estimated directly from the 


2 We note that MRLens offers a second method to determine 
significant wavelet coefficients b ased on the False Discovery Rate 
llStarck. Pires fe Refregiei] I 2 OO 6 I . and references therein), which 
may improve the quality of the reconstructions by limiting the 
ratio of false detections to true detections to a user-specified value. 
However, for ease of comparison with GLIMPSE, which defines 
significant coefficients based on their SNR, we chose to use the 
kcr method in this work. 
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noisy map, and we typically choose the threshold k to be 
between 3 and 5. Coefficients that are deemed thus to be 
significant are retained in the reconstruction, while an en¬ 
tropy based regularisation is applied to those coefficients 
that are not significant. An estimate of the convergence is 
built iteratively, through application of this regularisation 
on the residual between the data and the current conver¬ 
gence estimate. 


2.2 3D mapping with GLIMPSE 

The full 3D problem presents an additional difficulty com¬ 
pared to the 2D problem as we are trying to de-project the 
structures along the line of sight. However, the radial lensing 
operator Q is singular and leads to an ill-posed inverse prob¬ 
lem, which does not accept a single, stable, solution even in 
the absence of noise. A robust method to address this very 
general class of problems is sparse regularisation, as in the 
2D case, above. If the signal to recover is assumed to be 
sparse (only a small number of coefficients are non-zero) in 
an adapted dictionary, then a robust estimate of the original 
signal can be recovered by solving a optimisation problem 
with a penalty on the l\ no rm of the coefficients of the signa l 
in the chosen dictionary dStarck. Murtagh fc Fadilill20ld) . 
GLIMPSE implements this approach in the context of the 
3D mapping problem and recovers an estimate of the den¬ 
sity contrast S as the solution of the following minimisation 
problem: 

1 9 

min n II 7 - P 7K Q^a ||s +A || a ||i , (11) 

a Z 

where a are the coefficients of the estimate of the density 
contrast 5 in an appropriate dictionary <f> (i.e. S = 4?a), 
A is a parameter tuning the sparsity constraint and, as be¬ 
fore, E is the covariance matrix of the noise, assumed to 
be Gaussian. The first term in equation m is called the 
data fidelity term, and represents the £2 distance between 
the data and the current estimate of the solution, weighted 
by the noise covariance, while the second term represents a 
sparsity penalty, i.e. the equation above is minimised when 
the £\ norm of the vector of coefficients is as small as possible 
whilst being as close as possible to the input data. 

GLIMPSE uses a 2D-1D dictionary for $ composed of 
isotropic wavelets in the 2D angular domain and Dirac delta 
functions in the radial dimension. Such a representation is 
well adapted to the expected density contrast of a typical 
dark matter halo, which is always contained within a single 
redshift bin given the resolutions we are able to attain in 
the radial dimension. 

GLIMPSE solves this optimisation problem using a 
variant of the Fast Iterat i ve So ft Thresholding Algorithm 
fFISTA. iBeck fe Teboullel (120091 11. This algorithm relies on 
two steps, first the coefficients a of the solution are updated 
according to a gradient descent of the data fidelity term, 
then the coefficients are thresholded, setting to zero coeffi¬ 
cients of amplitude lower than A. This imposes the sparsity 
constraint on the solution. As with MRLens, GLIMPSE sets 
the threshold level to \j t n = kcrj tn where cry™ is the noise 
standard deviation at wavelet scale j and redshift bin n. 
With this choice of parameter, at each iteration only the 
significant coefficients (i.e. of amplitude above fcffj.n) are al¬ 
lowed to enter the solution. The only parameter to set for 



0.0 0.2 0.4 0.6 0.8 

Redshift 


Figure 1. This figure shows the virial mass and redshift for the 
96 clusters forming this study. 


the GLIMPSE algorithm is the significance threshold fc; as 
with MRLens, this is usually taken to be between 3 and 5. 


3 TESTS WITH CLUSTER SIMULATIONS 


In order to test the effectiveness of the two mapping tech¬ 
niques at detecting and reconstructing clusters of galax¬ 
ies, we simulated 96 cluster s of galaxies at a range 
of masses and redshifts as in iLeonard, Lanusse fe Starckl 
(120141) : Figure □ shows the points sampled in virial mass- 
redshift space. Each cluster ha lo follows an NFW profile 
dNavarro, Frenk fe Whitel 1 19971 ) with a concentration pa¬ 
rameter computed from the virial mass and redsh ift using 
the parameterisation given in lCoupon et al.l (|2012l ). 

From the NFW density profile, we computed the cor- 
responding shear sig nal, which is derived analytically in 
iTakada fc Jain! (120031 ). The computations required to simu¬ 
late the halos were performed making extensive use of the 
NICAEA software package^ using a flat ACDM cosmology 
with Om = 0.264, 9a = 0.736, and Ho = 71km/s/Mpc. 

Each cluster was simulated in an otherwise empty field 
with an angular pixel size of 1' x 1'. We considered our shear 
information to come from galaxies following a uniform an¬ 
gular distribution on the sky, with a redshift distribution 
given by 


n{z) = z a exp 



( 12 ) 


where we take zo = 1/1-4, giving a median reds hift z me d = 1, 
a = 2 and ft = 1.5 for a Euclid- l ike su rvey dTavlor et alJ 
120071 : iKitching. Heavens fc Miller! 1201 ll ). We assumed 30 
galaxies per square areminute, and that the intrinsic galaxy 
ellipticity follows a Gaussian distribution of width a e = 0.25. 


3 http://www2.iap.fr/users/kilbinge/nicaea/ 
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In their current forms, both MRLens and GLIMPSE take 
as inputs the shear, rather than the reduced shear (g = 
7/(1 — k)), and so our simulations produce noisy shear fields, 
rather than reduced shear fields. This is a simplification, in¬ 
cluded simply for ease of computation of our simulations, 
but which may lead to a biased recon struction if applied to 
real data. However, I Julio et al.l (I20l4i describe an iterative 
scheme to account for the reduced shear that converged af¬ 
ter three iterations and de-biases the reconstruction. Such 
an iterative process can be implemented with both MRLens 
and GLIMPSE in applications to real data to account for 
reduced-shear. This will be demonstrated for GLIMPSE in 
future work. 

Finally, we assumed Gaussian photometric redshift er¬ 
rors with a standard deviation given by a~ = 0.05(1 + z) 
and no bias. 

The 2D projected shear map was computed by inte¬ 
grating the lensing signal along each line of sight, while the 
3D shear map was computed using 30 tomographic bins of 
variable width A z, such that each tomographic bin con¬ 
tained the same number of galaxies. The bins were cho¬ 
sen to be sequential in photo-z space. Given the photomet¬ 
ric redshift errors, such a choice means that the true red- 
shi ft distributions for each bin ov erlap (see, e.g., figure 2 of 
iLeonard, Lanusse fe Starckll2014l 'l. Each map was generated 
with 64 x 64 angular pixels, each of width 56 = 1 arcmin. In 
order to properly asses the statistics of the reconstructions, 
for each cluster field, we computed 1000 noisy realisations 
of the data. 

3.1 Algorithm configurations 

In order to compare the two methods in a fair way, we used 
algorithm configurations that were as similar as possible. 
Both MRLens and G LIMPSE were run usi n g the sa me 2D 
wavelets (the starlet, IStarck. Fadili fc Murtaghll2007lj using 
6 wavelet scales in the decomposition. As GLIMPSE applies 
a ka hard threshold in density space, we configured MR¬ 
Lens to use ka thresholding, rather than the default False 
Discovery Rate (FDR) thresholding. This is not the optimal 
configuration for MRLens, as it may allow a higher frac¬ 
tion of false peaks to appear in a given reconstruction, but 
was chosen in order to provide a fair basis for comparison 
with GLIMPSE. The threshold level was chosen to be at 4<r 
for GLIMPSE. For MRLens, we considered three different 
thresholds: 3 ct, 3.5ct, and 4<r. 

In this paper, we are interested in both the probability 
of detecting a real cluster of a given mass and redshift with 
a given denoising threshold(the true detection rate), and the 
probability of finding a peak in a given field that rises above 
the denoising threshold but is due entirely to noise (the false 
detection rate). Comparing the true detection rate between 
2D reconstructions and 3D reconstructions at the same de¬ 
noising threshold will probe any intrinsic differences in the 
lensing signal to noise when the signal is projected in 2D vs 
when the full 3D information about the signal is retained. 
However, naively we expect the false detection rate at a fixed 
denoising threshold to scale with the number of pixels in the 
reconstruction, and therefore naturally expect a higher level 
of contamination from spurious peaks in a 3D reconstruc¬ 
tion of a given lensing field than in the 2D reconstruction of 
that field using the same denoising threshold. 


A mass reconstruction algorithm will ultimately be 
judged both on its ability to detect real clusters and its 
effectiveness at controlling spurious peaks. There is always 
a trade-off to be had between completeness, which may be 
increased by lowering the denoising threshold, and purity, 
which is improved by raising the denoising threshold. The 
setup of the present experiment does not allow us to com¬ 
pute the expected purity we might obtain in a reconstruction 
of a real weak lensing dataset. However, a comparison of the 
true detection rates obtained by MRLens and GLIMPSE 
when the denoising thresholds yield comparable false detec¬ 
tion rates provides a more fair assessment of the overall per¬ 
formance of the two methods in application to weak lensing 
studies. 

3.2 Peak identification 

We used CLFIND (IWilliams. de Geus fe Blitzlll994l 'l. which 
is a friends-of-friends algorithm, to identify connected pixels 
above a given threshold associated with a possible detection. 
The user is required to set a minimum value that a pixel 
can have in order to be considered as part of a detected 
structure. The aim is to identify the pixels associated with 
a detected peak out to large radii from the centre of that 
peak, so that we are able to accurately compute properties 
such as the density-weighted centroid for the detected peak. 

Setting the detection threshold in 3D is fairly straight¬ 
forward: the GLIMPSE algorithm applies thresholding to 
the solution, meaning that regions of the reconstruction in 
which there are deemed to be no significant coefficients dur¬ 
ing the iterative reconstruction process will be set to zero 
in the reconstruction. Therefore, any deviations from zero 
may be considered a detection of a feature above the noise. 
However, setting the CLFIND threshold to be zero can lead 
to blending of detected structures and an increased runtime 
for the algorithm, so to minimise any potential blending and 
to optimise the runtime of the CLFIND algorithm in 3D, we 
set the detection limit to be 5 m in = 1. 

Defining a detection in 2D is somewhat more compli¬ 
cated. While MRLens does apply hard thresholding during 
the iterative procedure, due to differences in the details of 
the algorithm - specifically, that MRLens computes the mul¬ 
tiresolution support and applies the entropy based regular- 
isation to the residual, which is then added back on each 
iteration to build the solution - there are no regions in the 
reconstructed convergence maps that are identically zero. 
Furthermore, due to the mass-sheet degeneracy, the overall 
normalisation of the reconstructed convergence map may be 
somewhat arbitrary. In practice, this is usually constrained 
by either setting the mean convergence (over a sufficiently 
large field) to be zero or, in cluster fields, by setting the 
convergence to be zero sufficiently far away from the clus¬ 
ter. We apply no such boundary conditions in the 2D re¬ 
constructions presented here. However, we note that in 3D 
the density contrast is naturally physically constrained to 
be ^ — 1 (i.e. the density is always positive), and GLIMPSE 
applies a positivity constraint to its estimate of the solution 
at each iteration. 

To ensure that we all the structures identified in the 
multiresolution support and no others, we perform a multi- 
step detection procedure on all the 2D reconstructions (see 
figure [2] for illustrations of each step): 
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(i) Identify all pixels in the convergence reconstruction 
(figure [2( a)) where the multiresolution support (figure [2(b) | ) 
is non-zero; find the minimum value of the convergence 
amongst these pixels, K m in■ We define Ki vm = min[3 x 

10 —3 , Kmin/2]. 

(ii) Search for all pixels with values Js Ku m that are con¬ 
nected to the pixels with non-zero multiresolution support. 

(iii) Mask all other pixels (figure 2(c) I. 

(iv) Run CLFIND on the masked map to identify and 
separate the peaks detected in the map (figure 2(d) I 


Note that the limiting value of 3 x 10 -3 above was chosen 
as an effective trade-off between capturing as many pixels 
as possible associated with each peak detected and minimis¬ 
ing blending and the CLFIND algorithm runtime. In cases 
where the convergence value at the locations of non-zero 
multiresolution support are close to this value, the limit is 
lowered in order to link a larger number of pixels to the de¬ 
tected peaks. This also ensures that peaks coinciding with 
non-zero multiresolution support are always detected. This 
does occasionally lead to blending, however CLFIND as the 
final step is able to isolate and separate blended peaks where 
two clear local maxima occur. 

We compute the location of each peak in 3D by comput¬ 
ing the density-weighted centroid [x,y,z], and in 2D anal¬ 
ogously by computing the convergence-weighted centroid 
[x,y\ using those pixels identified by CLFIND. 


3.3 Peaks arising from the noise 

Naively, we expect that when we consider the angular dis¬ 
tribution of all peaks detected in all 1000 reconstructions of 
a given cluster field, we will find a strong tendency to de¬ 
tect a peak in the centre of the field, where the cluster is 
located, and a roughly uniform distribution of false peaks 
arising from the noise, at a much lower amplitude. This is 
because false detections due to random noise are expected to 
be uniformly randomly distributed, as the noise is uniform 
and uncorrelated in the x — y plane. 

In figure [3] we plot a histogram of the x and y positions 
of all detected peaks in 1000 GLIMPSE 3D reconstructions 
of a field containing a relatively high signal-to-noise cluster 
of virial mass M v i r = 7 x 10 14 C 1 Mq located at a redshift 
of z c i = 0.35, and an analogous plot for reconstructions of 
a lower signal-to-noise cluster of M V i r = 9 x 10 13 /i _1 Mq 
at the same redshift. As expected, in both cases we see a 
concentration of detections around the centre of the image, 
and a roughly uniform background with the indication of 
some edge effects on the borders of the image. 

Figure [I] shows the comparable distribution of peaks 
from MRLens 2D reconstructions of the same fields. His¬ 
tograms are plotted for the 3 different denoising thresholds 
as indicated in the legend. In all cases, there is a concen¬ 
tration of detections around the centre, as expected, and a 
background of false peaks. The amplitude of both the cen¬ 
tral peak and the uniform background of false detections de¬ 
creases as the denoising threshold is raised. This is expected: 
the higher the threshold used, the more we will suppress the 
detection of peaks arising due to the noise. However, in do¬ 
ing this, we also reduce somewhat the probability that we 
will detect the central peak. 

When the edge effects at extremal values of x and y 



(a) The reconstructed con¬ 
vergence produced by MR¬ 
Lens, plotted on a logarithmic 
colour scale. 



(b) The multiresolution sup¬ 
port produced by MRLens, 
showing the pixel locations at 
which a significant wavelet co¬ 
efficient was detected. 



(c) Pixels not connected with 
regions of non-zero multireso¬ 
lution support are masked and 
CLFIND run on the result¬ 
ing image as illustrated above. 
We use here the same colour 
scheme as in figure[T(a) |above. 



(d) The output from 
CLFIND: a mask identi¬ 
fying and labelling distinct 
structures in the reconstruc¬ 
tion. The pixels associated 
with each clump are used 
to compute the convergence- 
weighted centroid of the 
peak. 


Figure 2. The pipeline for detection of peaks in 2D MRLens 
reconstructions . The 3D pipeline is simpler, involving application 
of CLFIND directly on the 3D reconstructions, with no masking 
required. 


are excluded, the distributions shown in figures [3] and 0] can 
be fit by a Gaussian plus a constant offset representing the 
uniform background of false detections. We carried out such 
fits for the distributions of x and y generated from the 1000 
reconstructions of each cluster field using GLIMPSE, and for 
MRLens for each of the three denoising thresholds tested. 
In all but the lowest signal-to-noise cases, where the fitting 
procedure failed, the mean x and y identified by the fitting 
procedure were clustered around \x,y\ = [32,32], and the 
standard deviations <j x and a y were consistently found to 
be in the range 0.3 — 1.4, with a trend to larger values of the 
standard deviation for fields with a lower signal to noise. 

Given the difficulty of obtaining a reliable Gaussian fit 
for cluster fields where the rate of detection of the cluster is 
very low, we chose to define a cluster detection in both 2D 
and 3D as any detection with 30 Sf x ^ 34 and 30 ^ y ^ 34 
in pixel units. Anything outside of this range is considered 
to be a false detection. This choice may artificially reduce 
the true detection count in lower signal-to-noise fields rel- 
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Figure 3. This figure shows the distribution of the x and y po¬ 
sitions of peaks detected in 1000 GLIMPSE 3D reconstructions 
of a field containing a cluster of virial mass 7 X 10 14 /i^ 1 Mq at 
a redshift of z c i = 0.35 (left panel) and the analogous plots for a 
field containing a smaller cluster of virial mass 9 X 10 13 /i — 1 Mq 
at the same redshift (right panel). In both cases, the central peak 
represents detections of the cluster itself, while the approximately 
uniform background represents false detections. Some edge effects 
are apparent at extreme values of x and y. 




Figure 4. This figure shows the distribution of the x positions 
(y positions are assumed to be symmetrical) of peaks detected in 
1000 MRLens 2D reconstructions of a field containing a cluster of 
virial mass 7 X 10 14 h — *Mg at a redshift of z c \ = 0.35 (left panel) 
and the analogous plots for a field containing a smaller cluster of 
virial mass 9 X 10 13 h —1 Mq at the same redshift (right panel). 
In both cases, the different coloured histograms are derived using 
different values of the denoising threshold. This choice clearly 
impacts both the false and true detection rates. 


ative to those with higher lensing signal-to-noise; however, 
in the majority of cases, we retain as true detections those 
peaks lying within 2 — 4 standard deviations from the mean 
position at the centre of the field. Moreover, the width of 
the x — y distributions for a given cluster field were found 
to be very similar regardless of the reconstruction method 
and denoising threshold. Therefore, the choice to restrict 
detections of the cluster to those falling within the range 
30 ^ x 34 & 30 ^ y ^ 34 does not unfairly (dis) advantage 
any particular method. 

We also note in the figures that there appears to be an 
overdensity of detections near the edges of the reconstruc¬ 
tions. In order to minimise the contamination from these 
edge effects in the analysis that follows, we exclude all false 
detections located within 4 pixels of the edge of the recon¬ 
struction in the transverse plane. 

In figure O we show a histogram of the mean false de¬ 
tection rate per reconstruction seen in each of our 96 cluster 
fields using GLIMPSE, and using MRLens with denoising 


Algorithm 

Threshold 

Mfit 


Median 

MAD 

MRLens 

3(7 

2.087 

0.357 

1.923 

0.231 

MRLens 

3.5cr 

0.373 

0.069 

0.377 

0.062 

MRLens 

4(7 

0.056 

0.018 

0.056 

0.018 

GLIMPSE 

4(7 

0.781 

0.047 

0.786 

0.044 


Table 1. This table gives the parameters of the best-fit Gaussian 
to the distribution of the number of false detections per realisation 
in the 96 cluster fields analysed with MRLens at 3 different de¬ 
noising thresholds and GLIMPSE at a denoising threshold of 4cr. 
Also listed are the median and median absolute deviation (MAD) 
statistics computed from the 96 cluster field false detection rates. 


thresholds of 3cr, 3.5cr and 4cr. Overplotted in the figures are 
the best-fit Gaussian distributions. The fit parameters, the 
mean fiat and the standard deviation aat are listed in Table 
|T]and compared with the median and median absolute de¬ 
viation (MAD) statistics computed from the 96 computed 
false detection rates. 

We can clearly see that the MRLens false detection rate 
in 2D is substantially smaller (by a factor of ~ 38) than that 
seen in 3D at the same denoising threshold of 4<r, and only 
becomes comparable to the 4<r GLIMPSE false detection 
rate when the denoising threshold drops between 3 — 3.5o\ 
As noted in section m given that there are more pixels in 
3D than in 2D, it is expected that we will have more false 
detections in 3D than in 2D for the same threshold level, so 
this result is unsurprising. 

On the other hand, it does indicate that when assess¬ 
ing the overall performance of the GLIMPSE and MRLens 
algorithms it is not sufficient to simply look at their clus¬ 
ter detection rates at the same denoising threshold. In both 
cases, there will be a trade-off between increasing the clus¬ 
ter detection rate by lowering the denoising threshold, and 
controlling the false detection rate with a stricter choice of 
threshold. However, in 2D the false detection rate is natu¬ 
rally lower than in 3D, and therefore a lower threshold can 
be chosen in 2D than in 3D for the same level of false peak 
contamination. 

For this reason, we compare the GLIMPSE results using 
a denoising threshold of 4<r with MRLens results at three 
different thresholds. Comparing the MRLens and GLIMPSE 
true detection rates when both use a 4cr threshold allows us 
to determine whether there is any intrinsic signal-to-noise 
advantage in studying lensing systems in 3D rather than 
2D. Any information gained in 3D comes at the price of 
increased false detections, and so we lower the threshold 
in 2D so that the false detection rates become comparable 
between GLIMPSE and MRLens, to determine whether the 
information gain in 3D is sufficient to outweigh the increased 
false detections compared with what is achievable in 2D. 


3.4 Control of false detections 

It is clear that spurious peaks arising from the noise may 
mimic real peaks in a reconstruction, and will therefore be a 
significant contaminant in any study that aims to constrain 
cosmological parameters using peak counts, or which aims 
to estimate the mass function. It is therefore important to 
consider strategies that might be implemented to remove 
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GLIMPSE false detection rate, 4a MRLens false detection rate, 3a 
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Figure 5. This figure shows the distribution of the mean false 
detection rate per reconstruction for the 96 cluster fields described 
in the text, reconstructed using GLIMPSE with a 4cr denoising 
threshold, and MRLens with denoising thresholds at 3cr, 3.5 cr, 
and 4 a. 



False detections per reconstruction 



False detections per reconstruction 
MRLens false detection rate, 3.5a 




these false detections or limit their impact on any resulting 
cosmological inference. 

Figures [3] and [4] indicate that the distribution of false 
peaks is uniform across the field under randomisations to 
the underlying noise. While this is true in our idealised sim¬ 
ulations, in which we have assumed a uniform distribution of 
background sources, real data often contains a complex gap 
structure and boundaries, which may give rise to a less uni¬ 
form distribution of false peaks. We nonetheless expect that 
peaks arising due to the noise will not appear consistently 
in the same places under different, randomised realisations 
of the noise. 

In practice, when applying either GLIMPSE or 
MRLens to real data, one would need to carry out 
Monte—Carlo or bootstrap resampling of the data in or¬ 
der to gain a full understanding of the noise properties 
of the reconstructed maps. More over, as demonstrated in 
iLeonard, Lanusse fc Starckl ll2014h . a reliable estimate of the 
mass and redshift of a detected cluster, and their associ¬ 
ated errors, can only be obtained by considering an ensem¬ 
ble of reconstructions carried out under randomisations to 
the noise in the data. 

Given that false detections are randomly distributed, 
the probability of detecting a noise peak repeatedly in the 
same position under a different realisation of the noise is 
very low. In other words, we expect the ’detection rate’ as¬ 
sociated with a given false detection at a position [*, y, z] 
to be very low. One way to visualise this is to consider a 
mean reconstruction, obtained by summing the reconstruc¬ 
tions obtained from the all 1000 realisations of the data. 

A real cluster should be detected consistently at the 
same location and with a similar peak amplitude in (many of 



Figure 6. The mean of 1000 reconstructions in 2D (3o- thresh¬ 
old; left panel) and 3D (4<r threshold; right panel) of field 53, 
containing a cluster of mass M v i r = 7 X 10 14 h ~ 1 M @ and red- 
shift z c i = 0.35. Both plots use a logarithmic colour scheme, to 
accentuate any low-level false detections in the mean reconstruc¬ 
tion. The detection rate for this cluster in both 2D and 3D at the 
chosen denoising thresholds is 100%. 



Figure 7. The mean of 1000 reconstructions in 2D (3.5 ct thresh¬ 
old; left panel) and 3D (4 ct threshold; right panel) of field 46, 
containing a cluster of mass M v i r = 9 X 10 13 /i _1 Mq and red- 
shift z c l = 0.35. Again, we use a logarithmic colour scheme. The 
detection rate for this cluster in both 2D and 3D at the chosen 
denoising thresholds is ~ 50%. 


the) reconstructions using different realisations of the noise, 
so this feature should be prominent in such a mean map, 
while the false detections, which appear in random positions, 
will be suppressed by averaging many reconstructions. 

Figure [6] shows the mean 2D (3.5 ct threshold) and 3D 
(4cr threshold) reconstructions of field 53. To highlight any 
low-level false detections, we plot the densities on a logarith¬ 
mic colour scale. While this particular cluster has a detection 
rate of ~ 100% in both 2D and 3D, the mean density maps 
are typical for clusters with a much lower detection rate (e.g. 
figure [7] where we show the results for field 46, containing a 
cluster of mass M v i r = 9 x 10 13 h~ 1 M@, redshift z c \ = 0.35, 
and detection rate ~ 25%). In 2D, it is trivial to identify the 
location of the true density peak in the mean reconstruction: 
no other significant peaks are visible in the mean map. 

In 3D, we notice several interesting features. Firstly, 
the cluster is clearly visible in the mean reconstruction, 
but appears to be smeared out in redshift. This is be¬ 
cause the broad lensing efficiency kernel Q gives rise to 
a large uncertainty in the redshift at which the cluster is 
detected. This is evidenced by the large error bars seen 
in the redshift estimates f or th ese clusters presented in 
ILeonard, Lanusse fe Starckl (120141) . The cluster is clearly de¬ 
fined, though, and no significant false detections are seen at 
redshifts below z < 1.5. 

We do see the significant presence of high-amplitude 
false detections at the liigh-redshift boundary of the 
mean reconstruct ion, h owever. Such fa lse det ections 
were seen in ILeonard, Dupe fe Starckl (l2012l 'l and 
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Redshift 

Figure 8. The redshift distribution of cluster detections and false 
peaks in reconstructions of 1000 noisy realisations of shear data 
for field 53, containing a cluster of mass M v i r = 7 X 1 0 1 4 h ~ 1 A'/q . 
redshift z c i = 0.35. 

iLeonard, Lanusse fc Starckl (l2014l h and arise due to 
overfitting of the shear data in the highest redshift bins. 

In figure 0 we plot the redshift distribution of the 
central cluster detections and the false detections for field 
53. While the redshift distribution of the cluster detections 
shows an approximately Gaussian distribution centred on 
the true input redshift (in this case, z c i = 0.35), the red¬ 
shift distribution of the false detections shows a significant 
overdensity at very low and very high redshift, and a flat 
distribution in between. This distribution is characteristic 
of all the cluster fields studied. Given that we physically 
do not expect to see such a high density peak at such high 
redshifts, truncating the 3D reconstructions at z = 1.5 will 
remove these high redshift false detections, which account 
for ~ 30% of all false detections seen in 3D, whilst having a 
minimal impact (< 2% on the cluster detection rate). 

3.5 Cluster detection rate: 2D vs 3D 

With false detections now well understood, we turn our 
attention to consider the detection rate expected for clus¬ 
ters as a function of mass and redshift for both MRLens 
and GLIMPSE. Figure [9] shows the fraction of noise realisa¬ 
tions in which a given cluster was detected as a function of 
mass for each of the 8 different redshifts at which we simu¬ 
lated cluster haloes. This figure effectively shows the selec¬ 
tion function for each method in the configurations specified 
above as a function of mass and redshift. Error bars were 
computed assuming Poisson noise only. 

The selection functions in 2D and 3D all show 
the same general trends: an asymptote towards zero 
detections at low mass and an asymptote at high 
mass towards 100% detections. There is also a con¬ 
sistent trend towards lower detection rates at high 
redshift than low redshift for fixed cluster mass. How- 


Weak lensing reconstructions in 2D & 3D 

ever, these dependencies differ between 2D and 3D: 

GLIMPSE (4(j) vs MRLens (4a): 

This comparison offers insight into whether there is a 
natural signal-to-noise advantage in studying weak lensing 
clusters in 3D rather than 2D. In 26 fields, both GLIMPSE 
and MRLens have detection rates of 100% (the GLIMPSE 
detection rate is 100% in 39 of the cluster fields). In every 
other field, GLIMPSE outperforms MRLens in terms of 
detection rate. The largest improvement is seen for field 88, 
a 2 x 10 14 /i _1 Mq cluster at a redshift of 0.75, where the 
GLIMPSE detection rate of 0.178 is a factor of 11.lx larger 
than that with MRLens (0.016). The GLIMPSE detection 
rate is 2x higher than that of MRLens in 27 fields and 
5x in 6 fields. There therefore appears to be a natural 
advantage to using 3D reconstructions for cluster detection 
in weak lensing surveys. 

However, as noted in section ED3 comparing the 
GLIMPSE and MRLens detection rates at the same de¬ 
tection threshold is not an entirely fair comparison, as the 
GLIMPSE reconstructions show a ~ 38 x higher rate of false 
detections than the MRLens 2D reconstructions at a 4a 
threshold. A more fair comparison can be obtained by reduc¬ 
ing the MRLens denoising threshold to a level such that the 
reconstructions show a similar rate of false detections. This 
level is reached at a threshold somewhere between 3 — 3.5a. 

GLIMPSE (4a) vs MRLens (3.5a): 

The GLIMPSE false detection rate at 4a is significantly 
closer (within a factor of ~ 2) to that of MRLens at 3.5a, so 
a comparison between these two sets of reconstructions is of 
particular interest. The detection rates in this case, partic¬ 
ularly at low redshift, are quite similar. In 31 cluster fields, 
both MRLens and GLIMPSE attain a 100% detection rate. 
In 17 cluster fields, the MRLens detection rate is higher 
than the GLIMPSE detection rate, but the difference is 
typically on the order of 5— 10%, with the largest difference, 
~ 20%, seen in field 1 where the GLIMPSE detection rate 
is 0.134 and the MRLens detection rate is 0.161. As the 
cluster is moved to higher redshift, the detection rates in 
3D become consistently higher than in 2D: in 48 fields, the 
GLIMPSE detection rate is higher than that of MRLens, 
and in 3 of those fields the improvement is greater than a 
factor of 2. Field 88 is again one of the fields in which the 
largest improvement is seen. In this case the GLIMPSE 
detection rate is ~ 2.7x larger than that of MRLens. 

GLIMPSE (4a) vs MRLens (3a): 

Given that the false detection rate of GLIMPSE at 4a 
remains significantly higher than that of MRLens at 3.5a, 
for completeness we again lower the MRLens denoising 
threshold to 3a. Here the MRLens false detection rate is 
~ 2.5x larger than that of GLIMPSE. With such a low 
denoising threshold in 2D, we see a corresponding boost 
in detection rates, as expected. In this case, we find that 
both GLIMPSE and MRLens attain a 100% detection rate 
in 34 cluster fields. In 44 of the cluster fields, the MRLens 
detection rate is higher than that of GLIMPSE. These are 
primarily at lower redshift, and the detection rates remain 
within a factor of 2 for all but 4 of the fields. The largest 
difference is seen in field 15, where the MRLens detection 
rate (0.224) is 2.7x larger than that of GLIMPSE (0.082). 
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Figure 9. The fraction of noise realisations in which the central cluster was detected as a function of mass and redshift. 3D mapping 
clearly outperforms 2D mapping for detection of clusters for all but the highest mass clusters at which both methods asymptote to a 
detection rate of 100%. 


At higher redshift, we find that the GLIMPSE detection 
rate remains higher than that of MRLens in 20 of the 
cluster fields, but the improvement here is moderate (the 
ratio between the GLIMPSE and MRLens detection rates 
is typically of order ~ 1.2). 


4 DISCUSSION 

The results above clearly demonstrate that mass mapping is 
a very useful tool for detecting clusters of galaxies in weak 


lensing surveys. At low redshift, GLIMPSE and MRLens 
provide reconstructions of comparable quality, when one 
considers both the cluster detection rate and the false detec¬ 
tion rate seen in these reconstructions. However, GLIMPSE 
mapping offers a distinct advantage in the detection of clus¬ 
ters at high redshift, showing a significantly higher sensitiv¬ 
ity to such clusters than seen in 2D with MRLens. 

The improvement seen in 3D is likely due to a combina¬ 
tion of effects. In projecting the data down into 2D, we lose 
some of the information content contained in the data. Ad¬ 
ditionally, particularly in the case of high redshift clusters, 
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the lensing signal becomes damped when averaging galaxy 
shapes over the full range of redshifts probed by the survey: 
in our experiment, the average cllipticity in each angular 
pixel [x, y] is given by 


e(as, 2 /) 


fo P{z)l(x,y,z)dz 
JgP{z)dz 


+ A/'(0,cr e /v / u) , 


(13) 


where 7 ( 2 , y, z) is in this case the analytically-derived shear 
at angular position [ x , y\ and redshift z, p(z) oc n(z) gives 
the probability distribution function for the galaxies in the 
survey and Af is a random noise contribution drawn from a 
Gaussian of mean 0 and standard deviation a e /y/n, where 
n is the number of galaxies per pixel. In practice, in a real 
galaxy survey, equation m will be replaced by a discrete 
sum over galaxy ellipticities. 

Clearly, the lensing signal 7 ( 2 ;, y, z) = 0 Mz < z c 1 , where 
Zd is the redshift of the cluster. If the cluster redshift is 
known a priori , one can retain only those galaxies behind the 
cluster; however in the absence of this knowledge, computing 
a 2D shear map by averaging over all the galaxies in a survey 
will result in a damping of the lensing signal due to the 
inclusion of galaxies that are not lensed, and this effect will 
be more pronounced for higher redshift clusters. 

The GLIMPSE 3D reconstruction method is signifi¬ 
cantly more sensitive than the state-of-the-art 2D recon¬ 
struction algorithm MRLens, particularly for clusters at low 
mass and/or at high redshift. At low redshift, the quality 
of GLIMPSE and MRLens reconstructions are comparable. 
The false detection rates are similar between GLIMPSE 
and MRLens at 3.5a, and the detection rates align simi¬ 
larly: GLIMPSE exhibits a detection fraction above > 50% 
for clusters of mass M v i r > 10 14 h 1 M @ at low redshift 
(z ~ 0.25); MRLens provides a comparable detection rate 
at a 3.5a at this redshift with the detection rates at 3a 
and 4a being 70% and 30%, respectively. At high red¬ 
shift, there is a clear improvement in 3D reconstructions. 
For example, GLIMPSE detects 63% of clusters of mass 
3 x 10 14 /i^ 1 Mq, while the MRLens detection rates are 
59% (3a), 37.5% (3.5a) and 19% (4a). At high redshift, 
where the weak lensing signal-to-noise is low, it is there¬ 
fore significantly advantageous to employ GLIMPSE for the 
detection of clusters in weak lensing maps. 

We have demonstrated that MRLens and GLIMPSE 
can be tuned to yield mass reconstructions with similar lev¬ 
els of purity; in other words, reconstructions in which the 
contamination due to peaks arising from the noise is at a 
similar level in both 2D and 3D reconstruction met hods. 
Given the findings of lPires. Leonard fe Starckl (|2012i 'l that 
peak counts undertaken on MRLens reconstructions provide 
a powerful method to probe non-Gaussianity and to discrim¬ 
inate between degenerate cosmological models, it is reason¬ 
able to expect that 3D reconstructions with GLIMPSE will 
offer an even more powerful constraint on non-Gaussianity, 
given the higher levels sensitivity to high-redshift clusters 
seen in GLIMPSE reconstructions. 

Furthermore, _as _ described in 

iLeonard, Lanusse fc Starckl d2014lh GLIMPSE recon¬ 
structions provide an unbiased estimator of both the 
masses and redshifts of the clusters detected. We expect 
that the selection function for optical halo finders will 
differ somewhat from that of GLIMPSE, but a detailed 
comparison would require a sophisticated simulation in¬ 


volving N-body cosmological simulations, raytracing, and 
a realistic semi-analytic modelling of galaxy formation. 
This comparison will be the subject of future work, and 
has implications for the study of galaxy clusters and their 
statistics in large upcoming surveys such as Euclid. 

A final cosmological application of GLIMPSE would 
be to measure the evolution of the high-mass end of the 
mass function. To do this, we need an accurate model of 
the GLIMPSE selection function as a function of cluster 
mass, redshift, and concentration parameter, and also as a 
function of the underlying cosmological model. Work is on¬ 
going in this area, and will be the subject of an upcoming 
publication. GLIMPSE is a non-parametric reconstruction 
method, which means that we do not compute a virial mass 
or redshift for the clusters we detect as part of the algo¬ 
rithm, but rather compute masses by integrating the density 
over a range of pixels in the reconstruction. In order to use 
GLIMPSE outputs to compute the cluster mass function, 
we must first be able to translate the GLIMPSE masses 
into virial masses, or to model the mass function for the 
GLIMPSE observables as a function of virial mass, redshift 
and cosmological model. This task is not trivial, but cer¬ 
tainly feasible. 

It is important to note that the simulations presented 
here are highly idealised, and that real data will pose many 
additional challenges. The data themselves will contain com¬ 
plicated shear and photometric redshift systematics, and 
the effects of more complicated (and possibly correlated) er¬ 
rors, biases, and photometric redshift catastrophic failures, 
as well as the impact of intrinsic alignments, on the resulting 
mass maps will need to be investigated. 

Moreover, a more realistic density distribution contain¬ 
ing clustered mass peaks will result in blending of struc¬ 
tures in mass maps, and thereby impact the measured peak 
counts and/or cluster mass fun c tion. In a very preliminary 
work, ILeonard, Dupe fc Stared (120121 ') demonstrated that 
sparsity-based 3D mapping techniques may be able to dis¬ 
entangle different density peaks located along the same line 
of sight, and while extensive tests have not yet been car¬ 
ried out with GLIMPSE, the indications are that GLIMPSE 
improves greatly on the results of this earlier work. A full 
exploration of this is ongoing, but this may highlight yet 
another advantage of using 3D reconstructions, rather than 
two-dimensional projections in which the line-of-sight infor¬ 
mation is lost. 

In addition, both magnification and instrumental effects 
will give rise to a spatially-varying survey depth, and in 
combination with the removal of foreground objects such 
as stars, this will give rise to a very complex survey mask 
that must be accounted for. It is clear that simple bin¬ 
averaging will not suffice in this case, and methods that can 
reconstruct the mas s distribution without resorting to bin¬ 
averaging (see, e.g. lLanuss 3 120151. , in prep) will be prefer¬ 
able. This represents a natural extension to the GLIMPSE 
approach, and this development is ongoing. 

Lastly, we point out that the experiment presented here 
was not blind, and the detection rates for clusters presented 
here do not represent the probability of detecting a cluster of 
a given mass and redshift in a weak lensign survey. Rather, 
they represent the probability of detecting a given cluster 
given that the cluster actually exists (i.e. they give the se¬ 
lection function for each reconstruction method). To extend 
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these results to predict the number of clusters we expect to 
detect in a given survey, one would need to convolve the 
cluster selection functions presented here with the cluster 
mass function, which gives the probability of the existence 
of a cluster of a given mass at a given redshift. 

Nonetheless, mass mapping remains in our view a po¬ 
tentially very useful tool for weak lensing cosmology, and 
GLIMPSE clearly offers several distinct advantages over 2D 
weak lensing reconstruction methods. Moreover, it provides 
a cluster detection method complementary to optical stud¬ 
ies. There is much work still to be done, but GLIMPSE 
is clearly a very promising tool for constraining the non- 
Gaussian part of cosmological density held and probing the 
evolution of nonlinear structures in the Universe. 


5 ACKNOWLEDGMENTS 

The authors would like to thank Filipe Abdalla and Paniez 
Paykari for useful discussions. This work is supported by 
the European Research Council grant SparseAstro (ERC- 
228261). We also gratefully acknowledge the thoughtful 
and insightful comments provided by the anonymous ref¬ 
eree, which have substantially improved the quality of this 
manuscript. 


REFERENCES 

Beck A., Teboulle M., 2009, SIAM J. Img. Sci., 2, 183 
Berge J., Amara A., Refregier A., 2010, ApJ, 712, 992 
Bernardeau F., van Waerbeke L., Mcllier Y., 1997, A&A, 
322, 1 

Borgani S., Guzzo L., 2001, Nature, 409, 39 
Coupon J. et al., 2012, A&A, 542, A5 
Dietrich J. P., Hartlap J., 2010, MNRAS, 402, 1049 
Falco M., Hansen S. H., Wojtak R., Brinckmann T., Lind- 
holmer M., Pandolb S., 2014, MNRAS, 442, 1887 
Girardi M., Giuricin G., Mardirossian F., Mezzetti M., 
Boschin W., 1998, ApJ, 505, 74 
Ivezic Z., Tyson J. A., Allsman R., Andrew J., Angel R.., 
for the LSST Collaboration, 2008, ArXiv: 0805.2366 
Jarvis M., Bernstein G., Jain B., 2004, MNRAS, 352, 338 
Julio E., Pires S., Jauzac M., Kneib J.-P., 2014, MNRAS, 
437, 3969 

Kaiser N., Squires G., 1993, The Astrophysical Journal, 
404, 441 

Kilbinger M., Schneider P., 2005, A&A, 442, 69 
Kitching T. D., Heavens A. F., Miller L., 2011, MNRAS, 
426 

Lanusse, F. e. a., 2015, in prep. 

Leonard A., Dupe F.-X., Starck J.-L., 2012, A&A, 539, A85 
Leonard A., Lanusse F., Starck J.-L., 2014, MNRAS, 440, 
1281 

Leonard A., Pires S., Starck J.-L., 2012, MNRAS, 423, 3405 
Lokas E. L., Mamon G. A., 2003, MNRAS, 343, 401 
Lokas E. L., Wojtak R., Gottlober S., Mamon G. A., Prada 
F., 2006, MNRAS, 367, 1463 

LSST Science Collaboration et al., 2009, ArXiv: 0912.0201 
Mamon G. A., Boue G., 2010, MNRAS, 401, 2433 
Moscardini L., Matarrese S., Mo H. J., 2001, MNRAS, 327, 
422 


Navarro J. F., Frenk C. S., White S. D. M., 1997, ApJ, 490, 
493 

Pires S., Leonard A., Starck J.-L., 2012, MNRAS, 423, 983 
Pires S., Starck J.-L., Amara A., Refregier A., Teyssier R., 

2009, A&A, 505, 969 

Refregier A., Amara A., Kitching T. D., Rassat A., 
Scaramella R., Weller J., Euclid Imaging Consortium f. t., 

2010, ArXiv: 1001.0061 

Rosati P., Borgani S., Norman C., 2002, ARA&A, 40, 539 
Schneider P., van Waerbeke L., Jain B., Kruse G., 1998, 
MNRAS, 296, 873 
Simon P., 2012, ArXiv/1203.6205 
Simon P. et al., 2011, MNRAS, 1789 
Simon P., Taylor A. N., Hartlap J., 2009, MNRAS, 399, 48 
Starck J.-L., Fadili J., Murtagh F., 2007, IEEE Transac¬ 
tions on Image Processing, 16, 2007 
Starck J.-L., Murtagh F., Fadili J. M., 2010, Sparse Image 
and Signal Processing. Cambridge University Press 
Starck J.-L., Pires S., Refregier A., 2006, A&A, 451, 1139 
Takada M., Jain B., 2003, MNRAS, 344, 857 
—, 2004, MNRAS, 348, 897 

Taylor A. N., Kitching T. D., Bacon D. J., Heavens A. F., 
2007, MNRAS, 374, 1377 

The Dark Energy Survey Collaboration, 2005, ArXiv As¬ 
trophysics e-prints 

VanderPlas J. T., Connolly A. J., Jain B., Jarvis M., 2011, 
ApJ, 727, 118 

Voit G. M., 2005, Reviews of Modern Physics, 77, 207 
Williams J. P., de Geus E. J., Blitz L., 1994, ApJ, 428, 693 


APPENDIX A: SUMMARY OF RESULTS 

The table below summarises the results of this study as pre¬ 
sented in §|33] Each simulated cluster held is denoted a held 
number, given in the hrst column of the table. Columns 2 
and 3 give the mass and redshift of the cluster in the held. 
The four left-most columns give the fraction of reconstruc¬ 
tions out of 1000 realisations of a noisy shear held in which 
the cluster was detected using MRLens with a denoising 
threshold of 3<r, 3.5<r and 4a, and using GLIMPSE with a 
denoising threshold of 4cr. 


Mass MRLens GLIMPSE 


Field 

(h~ 

'M 0 ) 

Z 

3cr 

3.5cr 

4cr 

4(7 

1 

3 

X 

10 13 

0.05 

0.338 

0.161 

0.051 

0.134 

2 

4 

X 

10 13 

0.05 

0.431 

0.252 

0.123 

0.221 

3 

5 

X 

10 13 

0.05 

0.579 

0.383 

0.214 

0.320 

4 

6 

X 

10 13 

0.05 

0.638 

0.447 

0.262 

0.417 

5 

7 

X 

10 13 

0.05 

0.738 

0.560 

0.350 

0.520 

6 

8 

X 

10 13 

0.05 

0.803 

0.664 

0.466 

0.610 

7 

9 

X 

10 13 

0.05 

0.830 

0.683 

0.507 

0.690 

8 

1 

X 

10 14 

0.05 

0.901 

0.795 

0.634 

0.756 

9 

2 

X 

10 14 

0.05 

0.991 

0.979 

0.949 

0.988 

10 

3 

X 

10 14 

0.05 

0.996 

0.997 

0.995 

1.000 

11 

4 

X 

10 14 

0.05 

0.999 

0.999 

0.999 

1.000 

12 

5 

X 

10 14 

0.05 

0.999 

0.998 

0.997 

1.000 

13 

6 

X 

10 14 

0.05 

1.000 

1.000 

1.000 

1.000 

14 

7 

X 

10 14 

0.05 

1.000 

1.000 

1.000 

1.000 

15 

3 

X 

10 13 

0.15 

0.224 

0.083 

0.024 

0.082 

16 

4 

X 

10 13 

0.15 

0.338 

0.156 

0.061 

0.149 


© 2014 RAS, MNRAS 000. ITIO 





Weak leasing reconstructions in 2D & 3D 13 


Mass MRLens GLIMPSE 


Field 

(h-'Me) 

z 

3cr 

3.5 a 

4(7 

4(7 

17 

5 

x 10 13 

0.15 

0.465 

0.265 

0.111 

0.231 

18 

6 

x 10 13 

0.15 

0.569 

0.351 

0.165 

0.335 

19 

7 

x 10 13 

0.15 

0.671 

0.459 

0.248 

0.442 

20 

8 

x 10 13 

0.15 

0.759 

0.555 

0.346 

0.545 

21 

9 

x 10 13 

0.15 

0.845 

0.658 

0.459 

0.641 

22 

1 

x 10 14 

0.15 

0.909 

0.762 

0.558 

0.732 

23 

2 

x 10 14 

0.15 

0.999 

0.996 

0.981 

0.996 

24 

3 

x 10 14 

0.15 

1.000 

1.000 

1.000 

1.000 

25 

4 

x 10 14 

0.15 

1.000 

1.000 

1.000 

1.000 

26 

5 

x 10 14 

0.15 

1.000 

1.000 

1.000 

1.000 

27 

6 

x 10 14 

0.15 

1.000 

1.000 

1.000 

1.000 

28 

7 

x 10 14 

0.15 

1.000 

1.000 

1.000 

1.000 

29 

6 

x 10 13 

0.25 

0.361 

0.198 

0.069 

0.192 

30 

7 

x 10 13 

0.25 

0.498 

0.270 

0.117 

0.267 

31 

8 

CO 

O 

X 

0.25 

0.586 

0.351 

0.158 

0.352 

32 

9 

x 10 13 

0.25 

0.636 

0.421 

0.220 

0.437 

33 

1 

x 10 14 

0.25 

0.707 

0.502 

0.290 

0.523 

34 

2 

x 10 14 

0.25 

0.998 

0.967 

0.905 

0.975 

35 

3 

x 10 14 

0.25 

1.000 

1.000 

0.995 

1.000 

36 

4 

x 10 14 

0.25 

1.000 

1.000 

1.000 

1.000 

37 

5 

x 10 14 

0.25 

1.000 

1.000 

1.000 

1.000 

38 

6 

x 10 14 

0.25 

1.000 

1.000 

1.000 

1.000 

39 

7 

x 10 14 

0.25 

1.000 

1.000 

1.000 

1.000 

40 

8 

x 10 14 

0.25 

1.000 

1.000 

1.000 

1.000 

41 

9 

x 10 14 

0.25 

1.000 

1.000 

1.000 

1.000 

42 

1 

x 10 15 

0.25 

1.000 

1.000 

1.000 

1.000 

43 

6 

CO 

O 

X 

0.35 

0.207 

0.098 

0.044 

0.110 

44 

7 

x 10 13 

0.35 

0.296 

0.129 

0.039 

0.155 

45 

8 

CO 

O 

X 

0.35 

0.373 

0.201 

0.080 

0.210 

46 

9 

CO 

O 

X 

0.35 

0.448 

0.254 

0.099 

0.267 

47 

1 

x 10 14 

0.35 

0.513 

0.289 

0.120 

0.333 

48 

2 

x 10 14 

0.35 

0.944 

0.831 

0.664 

0.887 

49 

3 

x 10 14 

0.35 

0.999 

0.988 

0.955 

0.995 

50 

4 

x 10 14 

0.35 

1.000 

1.000 

0.999 

1.000 

51 

5 

x 10 14 

0.35 

1.000 

1.000 

1.000 

1.000 

52 

6 

x 10 14 

0.35 

1.000 

1.000 

1.000 

1.000 

53 

7 

x 10 14 

0.35 

1.000 

1.000 

1.000 

1.000 

54 

8 

x 10 14 

0.35 

1.000 

1.000 

1.000 

1.000 

55 

9 

x 10 14 

0.35 

1.000 

1.000 

1.000 

1.000 

56 

1 

x 10 15 

0.35 

1.000 

1.000 

1.000 

1.000 

57 

9 

x 10 13 

0.45 

0.250 

0.117 

0.035 

0.157 

58 

1 

x 10 14 

0.45 

0.300 

0.137 

0.044 

0.199 

59 

2 

x 10 14 

0.45 

0.772 

0.574 

0.361 

0.710 

60 

3 

x 10 14 

0.45 

0.964 

0.889 

0.756 

0.959 

61 

4 

x 10 14 

0.45 

1.000 

0.988 

0.947 

0.997 

62 

5 

x 10 14 

0.45 

0.997 

0.997 

0.990 

1.000 

63 

6 

x 10 14 

0.45 

1.000 

1.000 

0.999 

1.000 

64 

7 

x 10 14 

0.45 

1.000 

1.000 

1.000 

1.000 

65 

8 

x 10 14 

0.45 

1.000 

1.000 

1.000 

1.000 

66 

9 

x 10 14 

0.45 

1.000 

1.000 

1.000 

1.000 

67 

1 

x 10 14 

0.55 

0.210 

0.094 

0.020 

0.115 

68 

2 

x 10 14 

0.55 

0.556 

0.322 

0.135 

0.491 

69 

3 

x 10 14 

0.55 

0.836 

0.681 

0.432 

0.834 

70 

4 

x 10 14 

0.55 

0.963 

0.896 

0.746 

0.969 

71 

5 

x 10 14 

0.55 

0.997 

0.975 

0.912 

0.997 

72 

6 

x 10 14 

0.55 

0.999 

0.991 

0.967 

1.000 

73 

7 

x 10 14 

0.55 

1.000 

1.000 

0.993 

1.000 

74 

8 

x 10 14 

0.55 

1.000 

1.000 

1.000 

1.000 

75 

9 

x 10 14 

0.55 

1.000 

1.000 

1.000 

1.000 

76 

1 

X 

o 

Oi 

0.55 

1.000 

1.000 

1.000 

1.000 


Mass MRLens GLIMPSE 


Field 

(h^Me) 

Z 

3(7 

3.5(7 

4(7 

4(7 

77 

1 

x 10 14 

0.65 

0.103 

0.037 

0.009 

0.067 

78 

2 

x 10 14 

0.65 

0.355 

0.186 

0.051 

0.309 

79 

3 

x 10 14 

0.65 

0.589 

0.375 

0.186 

0.634 

80 

4 

x 10 14 

0.65 

0.806 

0.620 

0.409 

0.862 

81 

5 

x 10 14 

0.65 

0.928 

0.806 

0.616 

0.963 

82 

6 

x 10 14 

0.65 

0.971 

0.921 

0.794 

0.994 

83 

7 

x 10 14 

0.65 

0.986 

0.968 

0.910 

0.999 

84 

8 

x 10 14 

0.65 

1.000 

0.988 

0.956 

1.000 

85 

9 

x 10 14 

0.65 

1.000 

1.000 

0.990 

1.000 

86 

1 

x 10 15 

0.65 

1.000 

0.998 

0.996 

1.000 

87 

1 

x 10 14 

0.75 

0.073 

0.015 

0.005 

0.041 

88 

2 

x 10 14 

0.75 

0.160 

0.065 

0.016 

0.178 

89 

3 

x 10 14 

0.75 

0.368 

0.189 

0.080 

0.415 

90 

4 

x 10 14 

0.75 

0.577 

0.344 

0.187 

0.657 

91 

5 

x 10 14 

0.75 

0.691 

0.497 

0.283 

0.838 

92 

6 

x 10 14 

0.75 

0.825 

0.651 

0.430 

0.937 

93 

7 

x 10 14 

0.75 

0.905 

0.783 

0.596 

0.979 

94 

8 

x 10 14 

0.75 

0.959 

0.888 

0.722 

0.995 

95 

9 

x 10 14 

0.75 

0.980 

0.934 

0.813 

0.999 

96 

1 

x 10 15 

0.75 

0.986 

0.962 

0.888 

1.000 
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